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We adopt the continuum limit of a linear, isotropic, homogeneous, transparent, dispersion- 
negligible dielectric of refractive index n and examine the consequences of the effective speed of 
light in a stationary dielectric, c/n, for D'Alembert's principle and the Lagrange equations. The 
principles of dynamics in the dielectric-filled space are then applied to the electromagnetic La- 
grangian and we derive the Hamiltonian and Hamilton's equations of motion for the fields and then 
macroscopic Maxwell equations. The resulting theory is mathematically equivalent to the familiar 
Maxwell-Heaviside formulation of continuum electrodynamics, but differs in the way that the ma- 
terial contribution is treated. A direct derivation of the symmetric total energy-momentum tensor 
from the field tensor for a dielectric is used to demonstrate the utility of the theory by resolving the 
Abraham-Minkowski electromagnetic momentum controversy. 
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The dynamical equations of electromagnetism are 
known collectively as the Maxwell equations. In a mate- 
rial that responds linearly to electric and magnetic fields, 
the macroscopic Maxwell-Heaviside equations take the 
form 
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in Heaviside-Lorentz units. Here, E is the electric field, 
B is the magnetic field, D = eE is the displacement field, 
H = B/ u is the auxiliary magnetic field, e is the electric 
permittivity, fj, is the magnetic permeability, and c is the 
vacuum speed of light. 

In this article, we derive equations of motion for macro- 
scopic fields from a Hamiltonian treatment of electrody- 
namics in a linear dielectric. We adopt the continuum 
limit of a stationary linear dielectric of refractive index 
n and examine the consequences of the effective light 
speed c/n on D'Alembert's principle and the Lagrange 
equations. We derive the Hamiltonian and Hamilton's 
equations of motion for electromagnetic fields in a di- 
electric from the Lagrangian density. From Hamilton's 
equations, we derive macroscopic equations of motion 
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for the electric and magnetic fields in a stationary dielec- 
tric medium, where B = V x A is the magnetic field and 
II = (n/c)dA/dt is the conjugate momentum field. With 
the exception of Poisson's Law, Eq. (I1.2dp . the equa- 
tions of motion for the macroscopic fields, Eqs. (11.21) . 
are mathematically equivalent to the Maxwell-Heaviside 
equations, Eqs. (|l.lj) . and Poisson's Law conforms in the 
limit of static fields or the Coulomb gauge. As the success 
of a new theoretical construct is often gauged by its abil- 
ity to resolve previously intractable problems, we use the 
variant form of Maxwell's equations, Eqs. (|1.2p . to derive 
the total electromagnetic energy-momentum tensor for 
the field in a dielectric, a result that has been sought for 
over a century. Thus, the Abraham-Minkowski momen- 
tum controversy [l|, 0] is resolved in favor of a traceless 
symmetric total energy-momentum tensor that is based 
on the Gordon [3| form gc = _(B x II) /c of the electro- 
magnetic momentum density 



II. PARTICLE DYNAMICS IN A DIELECTRIC 
FILLED SPACE 

In the continuum limit, a distribution of particles is 
regarded as a continuous medium and a property of the 
particles can be represented by a property density that is 
a continuous function of the spatial and temporal coor- 
dinates. In continuum electrodynamics, a linear medium 
is a region of space in which light travels from a source 
to an observer at a speed of c/n. We consider space to 
be entirely filled with a linear, isotropic, homogeneous, 
transparent, dielectric operating in a regime in which dis- 
persion, electrostriction, and magnetostriction are negli- 
gible and, for convenience, we apply the term simple lin- 
ear dielectric to this medium. Although dielectrics in the 
real world are much more complicated than this model of 
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a simple linear dielectric, theoretical models of physical 
phenomena are often based on reducing the complexity 
of the real world. Including too many details can hinder 
our understanding, especially if the theory becomes ma- 
terial specific, is based on phenomenological descriptions 
of effects, or becomes mathematically intractable. 

In the rest frame of the simple linear dielectric medium, 
the constant refractive index n is the only property of a 
dielectric that is significant to the current problem. Here, 
we derive the dynamics of discrete particles as though the 
particles travel, without impediment, through a region of 
space that has the property that the speed of light is c/n. 
The dynamics of discrete particles is then generalized in 
the usual manner to derive the dynamics of electromag- 
netic fields in a dielectric. 

Let the rest frame of the dielectric be S(t, x, y, z) with 
orthogonal axes x, y, and z. Then position vectors in S 
are denoted by x = (x,y,z). If a light pulse is emitted 
from the origin at time t — 0, then 



For a system of particles in equilibrium, the virtual 
work of the applied forces fi vanishes and the virtual 
work on each particle vanishes leading to the principle of 
virtual work 
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we can write D'Alembert's principle, Eq. (|2.9j) . as 
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describes wavefronts in the S system. Writing time as 

a spatial coordinate ct/n, the four-vector (ci/n, x) = in terms of the generalized forces 
(ct/n, x, y, z) represents the position of a point as a mat- 
ter of geometry Then, we need to adapt the usual 
theory of dynamics of particles in the vacuum [6j, [7[ to 
this coordinate system. 

For a system of particles, the transformation of the 
position vector Xj of the i th particle to J independent 
generalized coordinates is 
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If the generalized forces come from a generalized scalar 
potential function V @, then we can write Lagrange 
equations of motion 
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where t = t/n. Applying the chain rule, we obtain the 
virtual displacement 
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where L = T — V is the Lagrangian. The canonical mo- 
mentum is therefore 



and the velocity 
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of the i th particle in the new coordinate system. Substi- 
tution of 



in a linear medium. 

The field theory @, S] is based on a generalization of 
the discrete case in which the dynamics are derived from 
a Lagrangian density C. The generalization of the La- 
grange equation, Eq. (|2 . 1 3|) , for fields in a linear medium 
is 
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into the identity 
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where xq — ct/n is the time-like coordinate in the ma- 
terial and xi, X2, and X3 correspond to the respective x, 
y and z coordinates. We adopt the typical conventions 
that Roman indices run from one to three, Greek indices 
run from zero to three, and di represents the operator 
d/dxi. The conjugate momentum field 
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is used to construct the Hamiltonian density 
from which Hamilton's equations of motion 

dA v on 
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are derived. 
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III. EQUATIONS OF MOTION FOR FIELDS IN 
A DIELECTRIC 

We consider a stationary dielectric block illuminated 
at normal incidence from the vacuum by a plane 
quasimonochromatic electromagnetic pulse. The lin- 
ear, isotropic, homogeneous, transparent, and dispersion- 
negligible dielectric is draped with a gradient-index an- 
tireflection coating. The dielectric is mechanically rigid 
and spatial variation of the refractive index is sufficiently 
smooth that reflection, and the associated radiation pres- 
sure, can be neglected. Then the refractive index is a 
smoothly varying real and time- independent function of 
position in a finite region of space. 

We take the Lagrangian density of the electromagnetic 
field in the medium to be 
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in the absence of charges. Applying Eq. (|2.16p . the 
canonical momentum field 
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is substituted into Eq. (I2.17[) in order to construct the 
Hamiltonian density 
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Hamilton's equations of motion in the dielectric 
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are obtained from the Hamiltonian density, Eq. (13.31) , us- 
ing Eqs. (|2.18p . It is straightforward to reduce a second- 
order differential equation to first-order equations. To 
that end, we introduce a new field variable 



B = V x A. 
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Substituting the definition of B into Eq. (|3.4bl) produces 
the Maxwell- Ampere Law, 
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Thus, the second-order Hamilton's equation, Eq. (|3.4b[) 
has been written as two first-order equations, Eqs. (|3.5[) 
and (|3.6[) , as promised. The curl of the other Hamilton's 
equation, Eq. (|3.4a|i . is 
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known as Faraday's Law. The divergence of Eqs. (|3.5 
and (|3.6|) produces 



and Gauss's Law 
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The Maxwell-Heaviside equations, Eqs. (jl.ip . arc 
based on the results of experimental studies of electricity 
and magnetism that were performed in the 18th and 19th 
centuries that were codified into physical laws. Bearing 
the names of Faraday, Ampere, and Gauss, not all of 
the laws of electricity and magnetism survived in their 
original form. Based on purely theoretical arguments of 
self-consistency, Maxwell @ modified the Ampere law to 
include the effect of time-dependent electric fields. The 
changes to the Maxwell-Ampere and Faraday laws that 
are derived here are not nearly so profound. The two ver- 
sions of these laws are mathematically equivalent, that is, 
they can be mutually transformed into each other using 
vector identities that commute a scalar function of posi- 
tion with the curl and divergence operators. However, it 
has been shown that these types of transformations alter 
the way various physical processes are expressed [l^, [H| ■ 
There are also differences in the way some quantities 
are calculated. For example, derivation of the Fresnel 
relations requires conservation of energy in addition to 
Stokes's Law. In the next section, we show that the re- 
formulation of the Maxwell equations allows us to resolve 
the century-old Abraham-Minkowski controversy [l], Q 
for the electromagnetic momentum in a linear dielectric. 



IV. FIELD AND ENERGY-MOMENTUM 
TENSORS 

In the Maxwell-Heavisidc formulation of classical con- 
tinuum electrodynamics, there are two pairs of fields, 
{E,B} and {D,H}, and two field tensors. Here, there 
is a single pair of fields {II, B} and a single field tensor. 
The field tensor, 
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is obtained in the usual way from 
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for homogeneous materials. 

The reduction to a single field tensor and a single pair 
of fields results in a considerable simplification for certain 
derivations. For example, the energy-momentum tensor 
is defined in terms of the field tensor by 12] 
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is the four-divergence operator in the dielectric medium. 
Continuity equations and the energy-momentum tensor 
in an inhomogeneous dielectric are treated in Ref. fllij . 
Finally, integrating over all space a, the quantities 
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is the Maxwell stress tensor and g a @ is the diagonal met- 
ric tensor with non-zero elements g 00 = 1 and g n = — 1. 

The form of the energy-momentum tensor has been 
debated for over a century [§]. The best known candi- 
dates are the 1908 Minkowski tensor and the 1909 
Abraham [l6[ tensor. Unlike previous candidates, T a P 
satisfies necessary conditions for an energy-momentum 
tensor. The array is traceless and symmetric. Valid con- 
tinuity equations [TiJ 
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are derived from the four-divergence dp,T a P = of T a ° 
[H, where 
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arc temporally invariant if the refractive index is a slowly 
varying function of position. The conservation of the 
total energy, Eq. (I4.8a|) . and total momentum, Eq. (|4.8bl) . 
can be viewed as uniqueness conditions for the elements 
of the total energy momentum tensor that represent the 
total energy density and total momentum density of an 
isolated system [14j . 

V. SUMMARY 



We have recast classical continuum electrodynamics 
into a region of space in which the speed of light is c/n, 
instead of c. In principle, nothing is changed because the 
equations of motion can be transformed to the macro- 
scopic Maxwell-Heaviside equations using vector identi- 
ties: the two set of equations of motion are mathemat- 
ically equivalent in dispersion-negligible transparent di- 
electrics. In practice, the new formulation introduces a 
new perspective on classical continuum electrodynamics. 
As a demonstration, we presented a one-line derivation 
of the total energy-momentum tensor, a result that has 
been sought for decades. This derivation provides affir- 
mation for the total energy-momentum tensor that was 
recently obtained by construction from continuity equa- 
tions [4[ and from conservation properties |14l |. 
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